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The diffractive photoproduction of vector mesons is usually described considering the two-gluon 
(Pomeron) exchange, non-diagonal parton distributions and the contribution of the real part to the 
cross section. In this letter we analyze the diffractive photoproduction of the T at HERA using an 
alternative model, the Color Evaporation Model (CEM), where the cross section is simply determined 
by the boson-gluon cross section and an assumption for the production of the colorless state. We 
verify that, similarly to the J/ty case, the HERA data for this process can be well described by the 
CEM. Moreover, we propose the analyzes of the ratio R = ctt/pj/v to discriminate between the 
distinct approaches. 
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The successful operation of the HERA ep collider has opened a new era of experimental and theoretical investigation 
into diffractive vector-meson photo- and leptoproduction. On the experimental side, the HERA accelerator extends 
the accessible energy range by more than one order of magnitude over previous experiments. The HERA data show 
that the cross sections for exclusive vector-meson production rise strongly with energy when compared to fixed-target 
experiments, if a hard scale is present in the process. On the theoretical side, vector-meson production has proven 
to be a very interesting process in which to test the interplay between the perturbative and nonperturbative regimes 
of QCD (For a review, see for example While in inclusive process, as open heavy flavor production, the cross 
section is described in terms of a perturbative term associated to the cross section of the partonic subprocess, and 
a nonperturbative term represented by the parton distributions, an analogous factorization of hard and soft physics 
does not apply to quarkonium production rates ||, which constitute a small fraction of the total open flavor cross 
section. This makes quarkonium physics a challenging field . 

The current picture in the literature to describe the diffractive photo- and leptoproduction of vector-mesons assumes 
that the color singlet property of the meson is enforced at the perturbative level by the two-gluon (Pomeron) exchange 
In this approach the amplitude, in the target rest frame, is factorized as a sequence of events very well separated 
in time: (i) the photon fluctuates into a quark-antiquark pair; (ii) the qq pair scatters on the proton target, and 
(iii) the scattered qq pair turns into a vector meson. The interaction is mediated by the exchange of two gluons in 
a color-singlet state. Moreover, the two gluon exchange amplitude can be shown to be proportional to the gluon 
distribution xg(x,Q ), with x = (M 2 + Q 2 )/W 2 and Q rs \{M 2 + Q 2 ), where W is the jp center of mass energy 
and M is the invariant mass of the qq system. In the case of the production of a heavy meson, the presence of 
the heavy meson mass ensures that perturbative QCD can be applied even in the photoproduction limit. This 
approach has been improved, particularly in the role of the vector meson light cone wave function describing 
reasonably the HERA J/^f data. When applied to the recent HERA T data p,|J, it only reproduces the data if new 
effects, which significantly contribute, are considered ( a ) tne non-diagonal parton distributions, which probe 

new nonperturbative information about hadrons and are a generalization of the conventional parton distributions 
(For a review see Ref. ||); and (b) the real part of the scattering amplitude. In Ref. Q, a strong correlation between 
the mass of the diffractively produced state and the energy dependence of the total cross section was found, implying 
a distinct energy dependence for the T and J/^ photoproduction. One of the main motivations for the study of 
diffractive photoproduction of vector mesons in the Pomeron model is the possibility of obtain a sensitive probe of 
the behavior of the gluon distribution, due to the quadratic dependence on the gluon distribution |^,[^] of the cross 
section in this model. 

An alternative view of the diffractive photoproduction process was proposed recently in Ref. [0, where the J/^f 
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photoproduction was analyzed using the color evaporation model (CEM) p2| . In this case the color singlet is not 
enforced at the perturbative level and the cross section for the process is given essentially by the boson-gluon cross 
section plus an assumption for formation of the colorless meson. In the CEM for the diffractive photoproduction 
the cross section is linearly proportional to the gluon distribution. These authors have obtained a parameter free 
prediction which describes very well the HERA J/ty data. 

In this letter we extend the application of the CEM for the recent HERA T data and verify that, using the 
parameters determined in T hadroproduction, this model reasonably describes the experimental data with no need 
to introduce a colorless (Pomeron) exchange at the perturbative level, as well as, non-diagonal parton distributions 
or the contribution of the real part of the scattering amplitude. Moreover, we propose the analyzes of the ratio 
R = (Jr/cj/* to discriminate between the distinct approaches. 

Let us start from a brief review of the Color Evaporation Model (CEM). One of the main uncertainties in the 
quarkonium production is related to the transition from the colored state to the colorless meson. Initially, the qq pair 
will in general be in a color octet state. It subsequently neutralizes its color and binds into a physical resonance. 
Color neutralization occurs by interaction with the surrounding color field. If we enforce that a colorless object is 
already present at the perturbative level then a hard interaction with the surrounding color field is assumed, and a 
minimal of two gluons should be exchanged with the proton (Pomeron models). However, if the colorless object is 
produced at the nonperturbative level, then there is not a minimal restriction in the number of gluons exchanged with 
the proton. The CEM provides a simple and general phenomenological approach to color neutralization (see also Ref. 
|p3| ). In CEM, quarkonium production is treated identically to open heavy quark production with exception that in 
the case of quarkonium, the invariant mass of the heavy quark pair is restricted to be below the open meson threshold, 
which is twice the mass of the lowest meson mass that can be formed with the heavy quark. For bottomonium the 
upper limit on the bb mass is then Ittiq. Moreover, the CEM assumes that the quarkonium dynamics is identical 
to all bottomonium states, although the bb pairs are typically produced at short distances in different color, angular 
momentum and spin states. The hadronization of the bottomonium states from the bb pairs is nonperturbative, usually 
involving the emission of one or more soft gluons. Depending on the quantum numbers of the initial bb pair and the 
final state bottomonium, a different matrix element is needed for the production of the bottomonium state. The 
average of these nonperturbative matrix elements are combined into the universal factor F[nJ ], which is process- 
and kinematics- independent and describes the probability that the bb pair binds to form a quarkonium T(nJ PC ) of 
given spin J, parity P, and a charge conjugation C. Once F has been fixed for each state (T, T' or T") the model 
successfully predicts the energy and momentum dependence jl4],|l5| . 

Considering the elastic T photoproduction at HERA, the CEM predicts that the cross section is given by 

a[T(nJ PC )} = F[nJ PC ] <f[bb] , (1) 

where the large distance factors F can be written in terms of the probability 1/9 to have a color singlet pair after 
the soft interactions, times the fractions pi of total bottomonium carried by the different states, in a similar way as 
considered in [fl2"f , where the corresponding p factors for charmonium are claimed to be universal. The short distance 
contribution is 
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Here a[bb] is the spin- and color-averaged cross section for open heavy quark production, which describes the boson- 
gluon fusion process 73 — > bb; M.t is the invariant mass of the bb pair, mj is the bottom quark mass and 2ms is 
the BB threshold. The differential cross section for the boson-gluon subprocess is well known (see Eq. (4) in Ref. 
jnj for the expression at leading order). This will contribute to elastic photoproduction, since in LO all energy of 
the photon is transferred to the bb pair. Higher order corrections will be mostly important for the inclusive and the 
inelastic cross section. Thus, the soft interactions have a double effect - to eliminate the color of the cc pair, allowing 
quarkonium production, and to allow elastic production with a single perturbative gluon exchange, by having soft 
gluon exchanges taking place in the non-perturbative part of the process. 

In CEM, the effect of these soft interactions is implicit in the non-perturbative factors. An explicit modeling for the 
soft gluon exchanges is done in the SCI model p3| , which is based essentially in the same general idea. In that model, 
by requiring either rapidity gaps or leading proton in the final state, it was obtained explicitly a good description of 
several diffractive processes, including diffractive J/tp and open beauty production in the Tevatron |lq] . This shows 
that soft gluon exchanges really might play a role in diffraction. 

The CEM has been used to predict the production of T in hadron and nuclear processes |fb5f , with the parameter 
F[nJ p ] determined from the quarkonium production experimental data at fixed-target energies. A remarkable 
feature of the CEM is that this model also accounts (within K factors) for the quarkonium production at the Tevatron 
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fuj . As the fixed target T data have generally given the sum of T, T' and T" production, due to the low mass 
resolution to clearly separate the peaks, only the global factor F[T + T' + T"] = 0.044 has been determined. This 
is an effective value which reflects both direct production and chain decays of higher mass states. As discussed in 
Ref. O], isolation of the direct production cross section for each Tj requires the detection of the radiated photons 
associated with chain decays, which is not currently available but might be possible at the LHC. Considering some 
assumptions related to the decay chain the following values for F in direct Tj production have been estimated jl7) : 
F d [T] = 0.023, F d [T'} = 0.02, F d [T"} = 0.0074, where the d superscript indicates the F for direct production. In 
our formalism the F factors can be related to the p factors by the simple relation Fi — l/9p,, where the 1/9 is the 
probability to have the bb pair in a color singlet state after the soft gluon exchanges, and pi are the universal factors 
with give the fractions of the onium cross section carried by the different onium states. Their corresponding values 
are, from above, 

p d [T] = 0.207 , p d [T'} = 0.18 , p d [T"} = 0.066 . (3) 

Once the free parameters have been determined in the hadronic processes, we can use the CEM to predict the T 
photoproduction at HERA. This also follows the assumption used in || that the production ratios of T, T' and 
T" are the same as those measured in hadron-hadron collisions. A comment is in order here. In Ref. |17f | the T 
hadroproduction has been calculated considering the next-to-leading order contributions for the cross section, which 
implies that the factor F determined from data does not contain perturbative contributions beyond leading order 
and can be considered an universal factor which describes the probability for quarkonium production. Moreover, 
the calculations in Ref. |17| have used the MRS D~ parton densities as input, but as demonstrated in Jlqj , where 
an update from the analyzes of the T suppression for the CMS detector was made, the T data can be equally well 



described using the GRV 94 LO |19| parton densities. 

The cross section a[bb] is computed at leading order using the GRV 94 LO |ll| parton densities with m& = 4.75 
GeV, tub = 5.279 GeV and the renormalization and factorization scales set to p = (Mqq) 1 / 2 [jl2 llj, with Q = b. 



In Fig. [l] we show our predictions for T photoproduction at HERA energies. Both the total T + T' + T" and direct 
T production are presented, since the cross sections measured by the ZEUS and HI collaborations did not select the 
direct T state, but rather the data were integrated over an interval of the p + ' p~ mass which includes at least the 
T, T' and T" resonances. We verify that the current experimental data does not allow to discriminate the distinct 
contributions. We emphasize that our results are completely parameter free and that the CEM model reasonably 
describes the scarce experimental data. 

The simplicity of the CEM strongly contrasts with the number of assumptions necessary in the Pomeron models 
to describe the same set of data. Here we proposed a signature to discriminate between these models. As we have 
quoted above, the Pomeron models predict a stronger growth in energy for the diffractive T photoproduction than 
in J /\& photoproduction, due to the strong correlation between the mass of the diffractivcly produced state and the 
energy dependence of the cross section. In contrast, in the CEM model the growth of the cross section is directly 
determined by the gluon distribution xg(x,p), where p is the factorization scale. Therefore, the energy dependence 
of the ratio 

RcEM = , (4) 

can be used to discriminate between models. 

As an intermediate step to calculate the mentioned rate, we show in Fig. ^ our previous calculation of <7j/^,, 
contrasted with the more recent HERA data §]. We see that these data can be reasonable explained without tuning 
the parameters, which were taken as the previously used in Ref. Hj. This gives one more clue that the CEM can be 
used to explain clastic J/ip photoproduction. 

In Fig. H we present the energy dependence of the ratio Rcem = 0"t/<L£/* calculated using the CEM, where dx 
represents the direct T production and we have used the results from Ref. (TlJ and from above to calculate crj/$. Our 
results show that this ratio is almost constant in the kinematic range of HERA, in contrast to the Pomeron model 
B, where a steep rise of the ratio is predicted [cf. Ref. M, Rp om oc W}. This ratio was also obtained in 0, where 
the value predicted agrees with data within errors for the MRS gluon, but underestimates the result for the steeper 
GRV gluon. In our case, the ar/pj/y ratio agrees with the lower bound observed by ZEUS || for both MRS and 
GRV gluon parameterizations. 

The original motivation of the study of diffractive photoproduction was the possibility to extract the gluon distribu- 
tion inside the proton. However, the dependence on the gluon distribution is one of the major differences between the 
descriptions of the diffractive photoproduction using either the Pomeron model or the CEM. Whereas the Pomeron 
model has a quadratic dependence on xg, this dependence is linear in the CEM. Our result demonstrates that before 
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to extract the gluon distribution from HERA T and J/^> data, one should determine the correct description for this 
process, for example by measuring the energy behaviour of the or /cj/w ratio. 

The CEM describes a large range of data in hadro- and photoproduction, as shown in Refs. |l^,|ll|,^0| . Using this 
model, in this letter we obtain a parameter free description of the elastic T photoproduction at HERA energies. We 
verify that this simple model reasonably describes the experimental data, similarly to the Pomcron models. As a 
distinct feature between these models, the CEM predicts a softer energy dependence of the ratio between the T and 
J /ty cross sections. Of course, when more precise data become available, the discrimination between these models 
should be possible, which will allow us to conclude whether the CEM is only an alternative phenomenological model for 
the current energies or it contains some underlying non-perturbative dynamics which is important in both diffractive 
and non-diffractive quarkonium production. In the latter case, more theoretical studies are necessary to understand 
the soft interactions in the process of color neutralization. 
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FIG. 1. Predictions of the Color Evaporation Model for T photoproduction at HERA. Cross section as a function of the 
photon-proton center-of-mass energy W. The data points are the ZEUS and HI results for the direct production of the T state, 
which should be compared with the full line. The total T + T' + T" contribution is also shown (dashed line). 
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FIG. 2. Predictions of the Color Evaporation Model for J/<P photoproduction at HERA. Cross section as a function of the 
photon-proton center-of-mass energy W. The data points are the new HI results for J/* production. 
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FIG. 3. Predictions of the Color Evaporation Model for the energy dependence of the ratio Rcem = in photoproduction 
at HERA. Data from ZEUS. 
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